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The  generation  of  acoustic  waves  by  a  modulated  thermal  source  is  examined  for  anisotropic 
materials.  A  wave  equation  is  developed  to  include  the  effect  of  a  thermoacoustic  source  for  the 
anisotropic  case.  The  dependence  of  the  thermoacoustic  source  term  upon  the  material  elastic 
constants  is  identified  for  a  thermal  source  varying  in  one  dimension.  This  dependence  is 
examined  for  several  classes  of  crystalline  anisotropy  and  it  is  found  that  thermoacoustic 
generation  varies  with  crystallographic  orientation.  The  directions  of  maximum  and  minimum 
L,  thermoacoustic  generation  are  not  consistent  for  a  given  crystal  class  and  are  dependent  upon  the 

Is  specific  values  of  the  elastic  constants  for  the  material. 
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INTRODUCTION 

It  has  long  been  known  that  a  modulated  thermal  source 
may  be  used  to  create  acoustic  waves.  The  effect  was  discov¬ 
ered  in  the  1880s  by  Alexander  Graham  Bell,1  who  studied 
the  audible  sound  generated  by  illuminating  matter  with  a 
periodically  interrupted  light.  The  basic  principle  of  ther¬ 
moacoustic  generation  within  a  solid  involves  the  coupling 
of  energy  from  thermal  expansion  and  contraction  into  an 
acoustic  wave.7  The  problem  of  thermoacoustic  generation 
at  the  surface  of  an  isotropic  solid  was  studied  by  White,' 
who  related  the  production  of  elastic  waves  to  thermal  and 
elastic  properties  of  the  solid.  More  recently,  thermoacous¬ 
tic  generation  at  the  surface  of  a  solid  has  been  used  as  a 
means  of  spectroscopic  characterization4  and  as  a  source  for 
imaging.'  Thermoacoustic  generation  is  a  particularly  useful 
technique  to  study  the  thermal  and  elastic  properties  of  thin 
layers  of  material  because  thermoacoustic  waves  can  be  gen¬ 
erated  within  a  layer  which  is  thinner  than  the  generated 
acoustic  wavelength.  In  contrast,  a  layer  several  acoustic 
wavelengths  thick  is  required  to  interact  effectively  with  ex¬ 
ternally  generated  acoustic  waves. 

This  paper  assesses  the  effect  of  elastic  anisotropy  on  the 
generation  of  acoustic  waves  by  a  modulated  thermal  source. 
It  will  be  shown  that  the  solid  elastic  constants  enter  into  the 
efficiency  of  thermoacoustic  generation  and  this  results  in  an 
orientation  dependence  for  thermoacoustic  generation  with¬ 
in  anisotropic  solids. 


I.  THERMOACOUSTIC  WAVE  EQUATION 

The  first  step  in  this  analysis  is  to  develop  a  wave  equa¬ 
tion  for  thermoacoustic  generation  in  anisotropic  materials. 
The  presence  of  a  locally  nonuniform  temperature  distribu¬ 
tion  within  the  material  produces  an  additional  strain  with¬ 
out  a  corresponding  change  in  stress.'  The  thermal  strain  S0 
is  proportional  to  the  temperature  increase  ©  and  to  the  lin¬ 
ear  thermal  expansion  coefficient  a.  This  strain  is  purely 
dilatational  and  can  be  written  in  abbreviated  subscript  nota¬ 
tion7  as 


S&  =  a  0 


rY 

1 

1 

0 

0 
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(1) 


where  a  and  ©  are  taken  to  be  scalar  quantities.  The  strain 
due  to  acoustic  waves  in  the  material  Sa  can  be  obtained 
from  the  difference  between  the  total  strain  5  and  the  tem¬ 
perature  associated  strain  , 

Su=S-S,y,  (2) 

giving  the  strain  displacement  relation, 

S„=V,u-S0,  (3) 

where  u  is  the  displacement  and  V,  is  the  symmetric  part  of 
the  displacement  gradient.  Recall  the  acoustic  equation  of 
motion 


V  T  =  p 


d  ~u 

IP 


and  the  constitutive  equation 
T  =  c  :  S„, 


(4) 


(5) 


where  T is  the  stress,  p  is  the  density,  c  is  the  stiffness  matrix, 
and  the  double  dot  product  represents  matrix  multiplica¬ 
tion.  Substituting  Sa  into  the  constitutive  equation  gives 


T = c  S  —  a  c :© 


(6) 


0. 


This  is  the  stress-strain  relation  for  thermoacoustic  genera¬ 
tion  and,  if  it  is  simplified  to  the  case  of  isotropic  media,  it 
becomes 
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which  is  called  the  Duhamel-Neumann  law,8  where  B  is  the 
bulk  modulus  of  elasticity. 

Combining  Eqs.  ( 3 )— ( 5 )  gives  the  thermoacoustic 
wave  equation9 


V  -(c :  V,n) 


=  V-(c:5e). 


(8) 

_ I 


II.  THERMOACOUSTIC  SOURCE  TERM 


The  right-hand  side  of  Eq.  ( 8 )  can  be  viewed  as  a  source 
term  due  to  the  nonuniform  temperature  distribution.  Ex¬ 
panding  this  term  gives 


rn 


V  -(c  :  S@  )  =  aV  •  c  :  0 
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(9) 


For  the  most  general  anisotropic  case,  there  are  21  indepen¬ 
dent  constants  in  the  c  matrix  (e.g.,  triclinic  crystal  class) 
and  the  thermoacoustic  source  term  can  be  written  as  fol¬ 
lows: 


V  -(c  :  S0  )  =  (  (CU  +  C, 2  +  Cj3  )  +  (C16  +  C26  +  C36)  +  (C15  +  C25  +  C3S))a* 


+ 


+ 
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Examining  the  components  of  the  thermoacoustic  source 
term  indicates  that  thermoacoustic  generation  is  enhanced 
for  large  thermal  gradients,  large  coefficients  of  thermal  ex¬ 
pansion,  and  large  values  for  the  elastic  constants.  In  addi¬ 
tion,  anisotropic  materials  allow  for  generation  of  acoustic 
waves  in  directions  perpendicular  to  the  thermal  gradients. 
Under  the  assumption  of  propagation  in  an  isotropic  medi¬ 
um,  Eq.  (10)  becomes 

V  (c  :S0)  =  3BaVQ.  (11) 

With  the  additional  assumption  of  purely  dilatational  dis¬ 
placements,  the  thermoacoustic  wave  equation  ( 8 )  becomes 

n  2 

cnV2u -p^^lBaVQ,  (12) 

dt 2 

which  can  also  be  directly  obtained  from  the  Dunhamel- 
Neumann  law. 

In  order  to  predict  the  efficiency  of  thermoacoustic  gen¬ 
eration  in  anisotropic  materials,  the  thermoacoustic  wave 
equation  (8)  will  be  examined  for  the  simple  case  of  a  ther¬ 
mal  gradient  in  one  direction.  In  particular,  assume  that  har¬ 
monic  thermal  energy  is  absorbed  at  the  boundary  between 
two  semi-infinite  media  with  a  resulting  thermal  distribution 
which  varies  in  the  x  direction  only.  The  thermoacoustic 
wave  equation  in  the  x  direction  is  then 

^  u * _ 1  d  k,  _  (£n_+_gi2_+£u)  <7®  ..i, 

dx2  v2  dt2  c,,  dx 

where  v  =  Jc^Tp 's  acoustic  velocity.  The  effect  of  the 
elastic  constants  on  the  thermoacoustic  source  term  can  now 
be  isolated  to  the  factor 

=  (Cii  +  c,2  +  c,3)/C|,.  (14) 

Likewise,  for  thermoacoustic  generation  due  to  thermal  gra- 


+  c34)  +  —  (c13  +  c23  +  cn)jaz.  (10) 

I - 

dients  confined  to  the  y  or  indirections,  the  elastic  portion  of 

the  source  terms  becomes 

Py  —  (CI2  +  C22  +  C2l)/C2l<  (15) 

Pz  =  (C|3  +  c23  +  c33)/c33.  ( 16) 

III.  THERMOACOUSTIC  GENERATION  IN  CRYSTALS 

The  thermoacoustic  generation  factor/?  is  given  in  Ta¬ 
bles  I-V  for  various  examples  from  the  orthorhombic,  tetra¬ 
gonal,  trigonal,  hexagonal,  and  cubic  crystal  classes.  The 
values  of  the  elastic  constants  used  to  calculate  /?  were  ob¬ 
tained  from  Auld.'°  The  orthorhombic  crystals  (see  Table 
I)  have  generation  factors Px,Py,  and/?z  which  are  indepen¬ 
dent  because  of  the  independence  of  the  elastic  constants 

^11  2^C22T^CJ.l 

and 

Cl2^Cl.li^C2.V 

For  the  three  orthorhombic  crystals  in  Table  1,  no  simple 
pattern  is  observed  for  the  thermoacoustic  generation  fac¬ 
tors  Px,Py,  and pz.  The  relative  magnitude  of/? in  the  three 
calculated  directions  is  different  in  each  case:  For  barium 
sodium  niobate,  Pz>p„>Py\  for  iodic  acid,  P  x  >  P  z  >  /?  v; 


TABLE  I.  Thermoacoustic  generation  factor  for  orthorhombic  crystals. 


Material 

0. 

0, 

0; 

Barium  sodium  niobate" 

1.64 

1.63 

1.76 

Iodic  acid* 

1.78 

1.32 

1.39 

Rochelle  salt’ 

2.16 

1.90 

1.88 

‘Piezoelectric  crystal,  values  are  calculated  using  rf 
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TABLE  11.  Thermoacoustic  generation  factor  for  tetragonal  crystals. 


Material 

Tirana 

0. 

Ammonium  dihydrogen  phosphate  ( ADP  )* 

1.39 

2.18 

Barium  titanate* 

2.20 

2.83 

Calcium  molybdate 

1.77 

1.62 

Indium 

2.80 

2.84 

Lead  molybdate 

2.11 

1.97 

Potassium  dihydrogen  phosphate  (KDP)* 

1.90 

2.01 

Rutile 

2.16 

1.58 

Tellurium  dioxide 

2.31 

1.41 

*  Piezoelectric  crystal,  values  are  calculated  using  ch 


and  for  Rochelle  salt,  0  x  >0V  >0 Likewise,  the  range  of 
variation  of  0  with  direction  is  different  in  each  case:  for 
Barium  sodium  niobate,  7%;  for  iodic  acid,  26%;  and  for 
Rochelle  salt,  13%.  All  three  of  the  orthorhombic  crystals  in 
Table  I  are  piezoelectric  and  the  stiffening  due  to  this  effect 
has  been  ignored  in  the  calculation  of  0  by  using  cE,  the 
elastic  stiffness  matrix.  The  inclusion  of  piezoelectic  effects 
in  0  could  be  accomplished  by  use  of  cD,  the  elastic  stiffness 
at  zero  electric  displacement. 1 1 

For  tetragonal,  trigonal,  and  hexagonal  crystals,  the 
thermoacoustic  generation  factors  are  equal  in  the  x  and  y 
directions  (/?,„„,  =0X  =  0y )  due  to  the  equivalence  of  the 
elastic  constants 

C|.  =  c22 

and 

c:i  —  c  if 

For  the  tetragonal  crystals  in  Table  II,  no  clear  pattern 
emerges  for  the  relative  magnitudes  of  /?„,„,  and  0Z.  The 
examples  shown  are  equally  divided  between  cases  where 
/?lnin,  >Pi  and 0.  >/?lrans.  A  substantial  difference  between 
the  relative  magnitude  of /?„,„,  and/?,  is  observed  for  telluri¬ 
um  dioxide  (39%),  ADP  (36%),  and  rutile  (27%),  indi¬ 
cating  that  for  these  crystals  the  magnitude  of  thermoacous¬ 
tic  generation  will  vary  significantly  with  orientation.  Also 
note  the  large  absolute  values  of  0Z  for  barium  titanate 
(2.83)  and  indium  (2.84). 

For  the  trigonal  crystals  in  Table  III,  there  is  no  consis¬ 
tent  pattern  for  the  relative  magnitudes  of/?lran,  and/?, .  The 
examples  show  that  larger  thermoacoustic  generation  may 
be  observed  in  either  direction.  However,  note  the  consis¬ 
tently  small  values  of/?  for  the  trigonal  crystals  and  in  parti¬ 
cular  the  small  0  for  quartz  (1.22). 

For  the  hexagonal  crystals  in  Table  IV,  in  each  case 
0 iran*  >Pz<  indicating  that  thermoacoustic  generation  is 


TABLE  III.  Thermoacoustic  generation  factor  for  trigonal  crystal!,. 


Material 

^tnn» 

0, 

Lithium  niobate* 

1  63 

1.73 

Lithium  tantalate* 

1.55 

1  69 

Quartz* 

1.22 

1.22 

Sapphire 

I  55 

1  46 

*  Piezoelectric  crystal,  values  are  calculated  using  cf  . 
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TABLE  IV.  Thermoacoustic  generation  factor  for  hexagonal  crystals. 


Material 

^tratu 

0, 

Beryllium  oxide* 

1.47 

1.36 

Cadmium  selenide* 

2.14 

1.94 

Cadmium  sulfide* 

2.20 

2.09 

Titanium,  crystal 

1.99 

1.76 

Zinc  oxide* 

2.08 

2.00 

Zinc  sulfide* 

2.09 

1.97 

*  Piezoelectric  crystal,  values  are  calculated  using  c*. 


consistently  greater  in  the  transverse  directions  relative  to 
the  z  direction.  However,  note  that  the  list  in  Table  IV  is  only 
a  sampling  of  hexagonal  crystals  and  others  may  deviate 
from  this  pattern.  Note  that  in  the  x-y  plane  of  a  hexagonal 
crystal,  elastic  wave  propagation  is  isotropic  resulting  in 
thermoacoustic  generation  which  is  independent  of  direc¬ 
tion. 

For  cubic  crystals,  the  thermoacoustic  generation  factor 
0  is  equal  in  the  x,  y,  and  z  directions  and  is  designated  0IOO. 
The  factor  0  can  be  calculated  for  other  directions  by  appli¬ 
cation  of  a  coordinate  transformation  to  the  c  matrix.12  Us¬ 
ing  a  45°  rotation  about  the  z  axis,  the  effective  elastic  con¬ 
stants  for  propagation  in  the  [110]  direction  of  a  cubic 
crystal  are  given  as  follows: 

fn  =  Cc, ,  +c12  +  2c44)/2,  c[2  =  (c,,  +c,2  -  2c44)/2, 

(17) 

C!l  =  CI2> 

and 

0\m  =  ( c'u  +  ^12  +  c[3  )/c[| .  (18) 


TABLE  V.  Thermoacoustic  generation  factor  for  cubic  crystals. 


Material 

P  HIO 

P\  |l> 

1  1 

Aluminum,  crystal 

2  14 

2.04 

2.01 

Barium  fluoride 

1  91 

1.90 

1.89 

Bismuth  germanate* 

1.47 

1  48 

1  48 

Bismuth  germanium  oxide* 

1  48 

1.80 

1  95 

Diamond 

1  49 

1.35 

1.31 

Europium  iron  garnet 

1.85 

1  82 

1.81 

Gallium  arsenide* 

1.91 

1  55 

1  46 

Gallium  phosphide* 

I  89 

1  55 

1  46 

Germanium 

1.75 

1  45 

1.37 

Gold,  crystal 

2.69 

2.34 

2.25 

Indium  antimonidc* 

209 

1.71 

1.61 

Indium  arsenide* 

2.09 

1.67 

1.57 

Indium  phosphide* 

2  13 

1  73 

1.62 

Iron,  crystal 

2  19 

1.70 

1  58 

Lithium  fluoride 

1.76 

1.40 

1.31 

Magnesium  oxide 

1  61 

1.38 

1.31 

Nickel,  crystal 

2.28 

1.76 

164 

Silicon 

1  77 

1.51 

1.44 

Silver,  crystal 

2  50 

2.01 

1.89 

Sodium  fluoride 

1.50 

1.64 

1.69 

Strontium  titanate 

1.64 

1  56 

1.54 

Tungsten,  crystal 

1  79 

1.79 

1.79 

Yttrium  aluminum  garnet  ( Y  AG ) 

1  66 

1  65 

1  64 

Yttrium  gallium  garnet 

1.81 

1.7: 

1.74 

Yttrium  iron  garnet  ( YIG ) 

1.83 

1  84 

1  84 

*  Piezoelectric  crystal,  values  are  calculated  using  cf 
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By  application  of  an  additional  rotation  about  they'  axis,  the 
elastic  constants  for  propagation  in  the  [  1 1 1  ]  direction  are 
found  to  be 

C\ i  —  (cn  +  2C|2  +  4c44)/2, 

c"2  =  (cu+2cn-2c44)/2,  (19) 

C|)  =Cy. 2, 

and/?, , ,  can  be  calculated  from  the  c"  values  as  above.  Using 
these  expressions,  the  values  of/?l00,  P ,  10,  and  /?, , ,  are  given 
in  Table  V  for  various  cubic  crystals.  For  most  of  the  exam¬ 
ples  in  Table  V,  the  relative  magnitude  of  the  generation 
factor  is  P\oo>Puo>P\\\<  indicating  that  thermoacoustic 
generation  is  greatest  along  an  axis  [  100]  of  the  cubic  crys¬ 
tals.  However,  in  four  cases  (bismuth  germanate,  bismuth 
germanium  oxide,  sodium  fluoride,  and  yttrium  iron  gar¬ 
net),  the  preferred  direction  for  thermoacoustic  generation 
is  [111]. 

Experimental  investigation  of  the  thermoacoustic  effect 
in  crystalline  materials  has  been  reported  by  Tam  and 
Leung”  and  by  Vladimirtsev  et  al.u  Tam  and  Leung  used 
photoacoustically  generated  sound  to  measure  the  ultra¬ 
sonic  velocity  as  a  function  of  orientation  in  a  sample  of 
polycrystalline  aluminum.  Although  their  paper  does  not 
mention  variation  in  the  amplitude  of  photoacoustic  genera¬ 
tion  with  crystal  orientation,  their  experimental  apparatus  is 
well  suited  to  measurements  of  this  kind.  Vladimirtsev  et  al. 
produced  photoacoustic  generation  in  piezoelectric  crystals 
and  observed  the  excitation  of  normal  modes.  They  observed 
that  surface  optical  defects  increased  the  intensity  of  photo¬ 
acoustic  generation;  however,  they  do  not  mention  variation 
in  photoacoustic  generation  with  crystal  orientation. 

IV.  SUMMARY 

A  wave  equation  has  been  derived  which  includes  the 
generation  of  elastic  waves  by  a  modulated  thermal  source. 
The  dependence  of  the  thermoacoustic  generation  on  mate¬ 
rial  elastic  constants  has  been  explored  for  anisotropic  mate¬ 
rials.  Although  substantial  variation  in  thermoacoustic  gen¬ 
eration  is  predicted  for  changes  in  orientation,  the  directions 
of  maximum  and  minimum  generation  are  not  consistent  for 


a  single  crystal  class  and  instead  depend  upon  the  specific 
values  of  the  elastic  constants. 
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